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ABSTRACT 

The  interaction  of  an  oblique,  overdriven  detonation  wave  with  a  vorticity  disturbance  is 
investigated  by  a  direct  two-dimensional  numerical  simulation  using  a  multi-domain,  finite- 
difference  solution  of  the  compressible  Euler  equations.  The  results  are  compared  to  those 
of  linear  theory,  which  predict  that  the  effect  of  exothermicity  on  the  interaction  is  relatively 
small  except  possibly  near  a  critical  angle  where  linear  theory  no  longer  holds.  It  is  found 
that  the  steady-state  computational  results  whenever  obtained  in  this  study  agree  with 
the  results  of  linear  theory.  However,  for  cases  with  incident  angle  near  the  critical  angle, 
moderate  disturbance  amplitudes,  and/or  sudden  transient  encounter  with  a  disturbance, 
the  effects  of  exothermicity  is  more  pronounced  than  predicted  by  linear  theory.  Finally,  it 
is  found  that  linear  theory  correctly  determines  the  critical  angle. 


^Research  was  supported  by  the  National  Aeronautics  and  Space  Adiririistiation  under  NASA  Contract 
No.  NAS1-1860S  while  the  authors  were  in  residence  at  the  Institute  fo-.  Computer  Applications  in  Science 
and  Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton,  VA  23665. 
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I  INTRODUCTION 

The  passage  of  a  weak  shear  disturbance  through  a  reactive  shock  wave,  or  deto¬ 
nation,  was  examined  by  Jackson,  Kapila,  and  Hussaini^.  Supersonic  engines  based  on 
oblique,  overdriven,  reacting  shock  waves  have  been  proposed  as  a  possible  alternative 
to  the  SCRAMJET  for  high-speed  propulsion^.  It  is  still  a  mattter  of  exploration  as  to 
whether  or  not  such  waves  can  be  stabilized.  Of  particiilar  interest  in  this  investigation  was 
the  effect  of  heat  release  on  the  refraction  and  amplification  of  the  vorticity  disturbance 
and  the  simultaneous  generation  of  acoustic  and  entropy  signals  behind  the  overdriven 
detonation.  The  detonation  was  assumed  to  be  at  an  angle  to  the  base  flow,  and  the  nor¬ 
mal  Mach  number  of  the  gas  ahead  of  the  detonation  front  was  taken  to  be  greater  than 
the  Mach  number  of  a  Chapman- Jouget  wave.  The  vorticity  disturbance  was  assumed  to 
be  a  small  amphtude,  planar,  shear  wave  with  wave  vector  parallel  to  the  base  flow  (i.e. 
transverse  disturbances).  There  exists  a  critical  angle,  dependent  upon  the  exothermicity 
of  reaction  and  the  overdrive,  such  that  the  relative  velocity  of  the  base  flow  behind  the 
front  is  subsonic  for  9  <  0^  and  supersonic  for  9  >  9^.  In  the  former  case,  the  amplitude 
of  the  generated  acoustic  disturbance  is  exponentially  decaying  behind  the  detonation; 
in  the  latter  case,  the  amplitude  is  constant.  The  critical  angle  approaches  zero  as  the 
Chapman-Jouget  limit  is  approached.  It  was  found  that  the  vorticity  was  significantly 
amplified  by  the  exothermicity  and  that  the  generated  acoustic  response  is  most  affected 
by  exothermicity  near  the  critical  angle.  Furthermore,  the  manner  in  which  the  shape  and 
structure  of  the  detonation  are  altered  by  the  disturbance  was  also  investigated. 

The  analysis  was  accomplished  by  considering  the  detonation  wave,  which  consists 
of  a  lead  shock,  an  induction  zone  and  a  fire  zone,  as  a  discontinuity  on  the  length  scale 
of  the  disturbance.  The  discontinuity  separates  an  unburnt  mixture  of  reactants  from  a 
burnt  mixture  of  reaction  products.  After  superimposing  a  vorticity  disturbance  in  the 
form  of  a  planar  shear  wave,  the  generalized  Rankine-Hugoniot  conditions  were  linearized 
about  the  base  flow  providing  the  conditions  to  determine  the  amplitudes  and  angles  of 
the  transmitted  vorticity  wave  and  the  generated  acoustic  and  entropy  waves.  The  effect 
of  the  disturbance  on  the  structure  of  the  detonation  was  found  by  considering  the  limit 
of  large  activation  energy,  E.  The  thickness  of  the  induction  zone  (the  region  between  the 
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lead  shock  and  the  fire  zone)  is  a  measure  of  the  thickness  of  the  detonation.  The  lead 
shock  and  the  fire  zone  were  treated  as  discontinuities  on  a  length  scale  comparable  to  the 
induction  zone  thickness.  The  equations  governing  the  perturbations  within  the  induction 
zone  due  to  disturbances  ahead  of  the  lead  shock  were  derived,  and  the  fire  zone  position 
was  determined  as  the  point  at  which  the  the  solutions  become  singular.  It  was  found  that 
disturbance  amplitudes  0{E~^)  have  an  0(1)  effect  on  the  fire  zone  position. 

The  aim  of  the  present  investigation  is  to  examine  the  predictions  of  linear  theory 
and  to  determine  by  means  of  a  direct  two-dimensional  numerical  simulation  the  regions  of 
validity  for  linear  theory.  For  this  investigation,  the  detonation  is  treated  as  a  discontinuity, 
and  the  generalized  Rankine-Hugoniot  relations  are  used  to  provide  appropriate  jump 
conditions  across  the  detonation.  The  effect  of  the  disturbances  on  the  internal  structure 
of  the  detonation  is  not  considered.  Of  particular  interest  is  the  behavior  of  the  solution 
for  angles  near  critical  where  linear  theory,  if  it  applies,  predicts  that  exothermicity  can 
have  a  significant  effect  on  the  generation  of  acoustic  signals. 


II  NUMERICAL  SIMULATION 

To  study  the  interaction  of  a  detonation  wave  with  a  vorticity  disturbance,  a  numerical 
approach  similar  to  the  one  utilized  by  Zang,  Hussaini,  and  Bushnell^'^  is  considered. 
The  major  differences  are  that  the  present  routine  takes  advantage  of  the  periodic  nature 
of  the  solution  in  the  transverse  direction  and  also  uses  a  multi-domain  scheme  in  the 
normal  direction.  Rather  than  using  coordinates  in  which  the  position  of  the  detonation 
is  stationary,  it  is  assumed  that  at  time  t  =  0  an  infinite,  planar  detonation  wave  starts  at 
X  =  0  and  propagates  into  the  unburnt  mixture.  The  position  of  the  detonation  x,  (y,  t)  is 
calculated  by  using  the  improved  shock-fitting  approach  of  Kopriva,  Zang  and  Hussaini^, 
generalized  to  include  exothermicity.  The  initial  conditions  are  chosen  such  that  in  the 
absence  of  disturbances  the  detonation  wave  will  propagate  to  the  right  with  a  Mach 
number  (relative  to  the  gas  in  front  of  the  detonation  wave)  greater  than  the  Chapman- 
Jouget  Mach  number.  The  base  flow  in  the  region  ahead  of  the  detonation  has  only  a 
vertical  component  so  that  the  net  result  of  the  velocity  of  the  detonation  front  in  the 
positive  x-direction  and  the  vertical  component  produce  a  base  flow  at  an  angle  d  in  a 


coordinate  system  attached  to  the  detonation  front.  Superimposed  on  the  uniform  base 
flow  is  a  vorticity  disturbance  which  propagates  at  the  same  angle  6  to  the  x-axis  (see  Figure 
1).  The  disturbance  is  assumed  to  have  unit  wavelength,  2n/\k\  =  1  where  k  is  the  wave 
vector.  The  flow  is  therefore  periodic  in  the  y-direction  with  period  27r/fcy  =  sec  6.  The 
initial  flow  behind  the  detonation  front  has  a  vertical  component  such  that  the  tangential 
velocity  across  the  front  is  continuous.  It  should  be  noted  that  the  flow  ahead  of  the  front  is 
prescribed  and  is  used  to  impose  the  appropriate  jump  conditions  across  the  discontinuity. 

The  physical  domain  in  which  the  fluid  motion  is  computed  is  given  by 

xi  <  X  <  x,{y,t),  0  <  y  <  sec(^),  t  >  0,  (1) 


where  x,  is  the  shock  position  and  x,{y,0)  =  0.  The  left  boundary  is  some  suitably 
chosen  negative  number  (usually  minus  one).  This  domain  takes  advantage  of  the  periodic 
nature  of  the  solution  in  the  y-direction.  In  the  x-direction,  a  multi-domain  approach  is 
used  which  allows  for  a  greater  number  of  grid  points  near  the  detonation  front  where 
greater  accuracy  is  needed  and  fewer  grid  points  away  from  the  front  where  the  solution 
is  smoother.  The  interfaces  between  the  domains  are  denoted  by  where  the  left 

boundary  is  xo{y,t)  =  Xi,  the  detonation  front  is  x„(y,t)  =  x,  (y,#),  and  the  planar 
interfaces  between  the  domains  art  Xi{y,t)  =  Xj(t),for/  =  l,n  — 1.  Since  the  behavior  of  the 
solution  near  the  critical  angle  is  of  interest,  it  is  necessary  to  continue  the  calculations  until 
the  detonation  wave  penetrates  a  significant  distance  into  the  unburnt  mixture.  Because 
of  this  expanding  domain,  a  greater  number  of  grid  points  in  the  x-direction  is  needed 
at  the  end  of  the  calculation  than  is  needed  in  the  beginning.  The  resolution  in  each 
subdomain  is  monitored,  and  the  number  of  grid  points  is  increased  when  necessary.  A 
linear  interpolation  is  used  to  transfer  the  solution  to  the  finer  grid.  The  numerical  routine 
also  allows  for  ramping  of  the  amplitude  of  the  imposed  disturbance,  thereby  decreasing 
in  size  the  portion  of  the  transient  caused  by  the  sudden  encounter  of  the  undisturbed 
detonation  with  the  vorticity  wave.  The  approach  to  steady  state  was  smoother,  allowing 
an  accurate  comparison  to  linear  theory. 

The  following  change  of  variables  is  made  to  computational  coordinates 

X  -  a:i_i(y,t) 


X,  = 


Xi{y,t)  -  Xi_i{y,t) 


i  =  l,n 


(2.a) 
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Y  =  y  cos(fl) 
T  =  t 


(2.b) 

(2.C) 


The  computationaJ  domains  are  therefore 

0  <  <  1,  0  <  y  <  1,  r  >  0. 


(3) 


The  fluid  motion  is  governed  by  the  two-dimensional  Euler  equations.  In  terms  of  the 
computational  coordinates  these  are 


where 


Qr  4-  BQxi  "I"  CQy  —  0  j  —  1,  n 


Q=[P,u,v,S\'^ 

Ui  7-Y.x 


B  = 


lu 

0 

0  0 


0 


c  = 


0 

0 


Cr  0 

0  Ui 


V 

7I; 

0  ■ 

c^yji 

V 

0 

0 

cXh 

0 

V 

0 

0 

0 

0 

V- 

i  «  +  uXi 

X  +  c 

A'i, 

y  =  y,  -I-  iiY,  -t-  vYy 


(4) 


(5a) 


(56) 


(5c) 

{5d) 

(5e) 


Here,  P  is  the  natural  logarithm  of  pressure,  c  is  the  local  sound  speed,  and  5  is  the 
entropy  divided  by  specific  heat  at  constant  volume.  The  velocities  u  and  v,  in  the  i  and 
y  directions  respectively,  are  scaled  by  the  sound  speed  ahead  of  the  shock.  The  ratio  of 
specific  heats  7  is  taken  to  be  1.4  for  all  calculations. 

As  in  Zang,  Hussaini  and  Bushnell3’4,  the  equations  are  discretized  in  each  subdomain 
using  the  finite  difference  method  of  MacCormack®,  which  is  a  variant  of  the  Lax-Wendroff 
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method.  Conditions  at  the  right  boundary,  at  the  interface  between  the  domains,  and  at 
the  left  boundary  need  to  be  imposed  to  determine  the  system  completely  .  At  the  right 
boundary,  the  improved  shock-fitting  routine  of  Kopriva,  Zang  and  Hussaini^  has  been 
generalized  to  include  exothermicity  and  time  dependence  of  the  flow  ahead  of  the  shock. 
Details  are  given  in  the  Appendix.  At  the  interface  between  domains,  a  routine  similar  to 
Kopriva®  is  used.  Since  the  interfaces  are  perpendicular  to  the  X-axis,  the  derivatives  in  Y 
can  be  calculated  by  using  the  values  of  Q  along  the  interface.  The  values  of  the  derivatives 
in  X,  however,  must  be  taken  from  the  left  or  right  of  the  interface  or  a  combination  of 
both.  The  method  chosen  here  is 

Ql  +  -  \BnQ^  +  =  0  (6) 

where  |51*  is  an  approximation  to  |5|  =  P|A|P"^,  A  is  a  diagonal  matrix  with  the  eigen¬ 
values  of  P  m  the  diagonal  and  P  is  a  matrix  of  the  eigenvectors  of  B.  The  superscript  I 
refers  to  values  calculated  on  the  interface  and  the  superscripts  R  and  L  refer  to  the  finite 
difference  approximations  to  the  X  derivative  in  the  subdomain  to  the  right  and  to  the 
left  of  the  interface  respectively.  The  approximation  used  fc;  rh^se  calculations  is 


\\U  +  C<t>\+\\U  - 

C(p\ 

0 

1 

1 

-1- 

0 

0 

\U\ 

0 

0 

\\U  +  c<t>\  -  - 

c<f>| 

0 

^\(j  -t-  c(p\  -l-  ^\U  —  C(p\ 

0 

0 

0 

0 

where  4>  =  (Xf The  conditions  imposed  at  the  left  boundary  naturally  depend  on 
the  inflow  at  the  boundary  being  supersonic  or  subsonic.  If  the  inflow  is  supersonic,  then 
all  four  components  of  Q  can  be  prescribed;  however,  if  the  inflow  is  subsonic,  only  three 
of  the  four  components  can  be  prescribed  at  the  left  boundary.  For  the  parameters  of  this 
study,  the  inflow  is  always  supersonic.  It  should  be  mentioned  that  the  normal  component 
of  the  flow  behind  the  detonation  is  supersonic  in  the  laboratory  frame;  however,  the 
normal  component  of  the  flow  behind  the  detonation  is  subsonic  in  a  frame  of  reference 
moving  with  the  detonation. 
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Ill  RESULTS  AND  COMPARISON  WITH  LINEAR  THEORY 


Lineal  Theory. 

Consider  a  coordinate  system  attached  to  the  detonation  front  such  that  the  speed  of 
the  detonation  is  zero.  The  flow  ahead  of  the  front  is  assumed  to  be  at  an  angle  6  to  the 
front,  the  normal  component  of  the  velocity  has  Mach  number  Mq,  and  the  vertical  com¬ 
ponent  of  the  velocity  is  VTq  (see  Figure  1).  The  generalized  Rankine-Hugoniot  conditions 
for  a  detonation  with  heat  release  parameter  a  provides  that  the  normal  component  of  the 
flow  behind  the  front  has  Mach  number  (relative  to  the  upstream  speed  of  sound) 


1 


rno 


1  +  iMq _ 

(7  +  l)Mo  (7  +  1) 


.u,?  “  ^1’  + 


1/2 


and  requires  that  the  tangential  components  be  continuous 

u’o  =  IVo- 


(8) 


(9) 


Therefore,  the  angle  of  the  base  flow  behind  the  detonation  is  given  by 

tan  <7  =  — ^  tan  (10) 

rno 

Superimposed  on  the  upstream  base  flow  is  a  sinusoidal  disturbance  in  the  velocity  given 
i>y 

=  (Uocos(k'  ■  x)  (11) 


where  Cq  =  (-Mq,  ITo)  and  Cq  ■  k'  —  0.  When  this  constant  pattern  of  vorticity  is  con- 
vected  through  the  detonation  front,  the  vorticity  w'ave  is  refracted  and  amplified,  and  an 
acoustic  and  entropy  wave  are  generated.  Linear  analysis  provides  analytic  expressions  for 
the  amplitudes  and  angles  of  all  three  disturbances  downstream.  The  result  is  that  the 
following  planar  waves  are  superimposed  on  the  base  flow  and  the  pressure  downstream 
{x  <  0): 


U1  =  1 

Ul 

0 

II 

cos(fc  •  x  +  6,) 

+  (\Uo\P{x) 

COs(fcp  •  X  -f  6p) 

.  VI  . 

1  1 

0 

■  (3  sin(fcp  ■  I  -f  (5p)  ■ 

Pi  =  -e7|17o|Afosec<^P(x)cos(fcp  ■  x  +  6^) 


(13) 
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where 


f  50  M  >  1 

5=1 

5o  n  <  1, 


^=1 

(,  rbexp(Ax)  M  <  1, 


(14) 


(15) 


and  jj.  is  the  local  Mach  number  of  the  relative  base  flow  behind  the  detonation  front.  The 
constants  S®,  So,  P®,  Po,  and  A  are  given  in  Jackson,  Kapila,  and  Hussaini^.  The 

vector  k  is  perpendicular  to  the  downstream  base  flow  and  the  vector  is  perpendicular  to 
the  direction  of  acoustic  propagation  which  makes  an  angle  (p'  with  the  x-axis  (see  Figure 
2).  This  angle  is  given  by 


r 


1.  —  tan 


0  —  cot  ^(,S) 
[(///dn)^cos^0tano] 


p  >  1 
H  <  1 


(16) 


with  3  =  (|1  -  and  =  1  —  j.c'^  cos^  0.  A  similar  disturbance  is  also  superimposed 

on  the  temperature  and  density  due  to  the  generated  entropy  wave. 


Nonlinear  Calculations. 

The  comparison  of  the  two-dimensional  calculations  to  the  linear  theory  consists  of 
comparing  the  amplitudes  of  the  vorticity  and  acoustic  waves  behind  the  detonations  with 
the  predictions.  To  achieve  this  comparison,  a  least  squares  fit  to  the  functional  forms 


Uy  -  u,  =  Sicos(A;'y)  -|-  52sin(A:^y)  (17) 

P  -  Po+  Picos(fc'y)  +  P2sin(^'y)  (18) 

is  performed  for  the  computed  vorticity  and  acoustic  fields  for  each  value  of  x.  The 
amplitudes  of  the  above  vorticity  and  acoustic  disturbances  5  =  (S^  -b  and  P  = 

(Pi  4-  are  divided  by  €\Uo\k'^  and  e7|^o|,W’osec<^  respectively  and  graphed  as  a 

function  of  x  for  various  times  t.  In  the  calculations  presented  here,  the  front  is  initially 
at  position  x  =  0,  and  behind  the  front  there  is  no  disturbance  (i.e.  P  =  0  and  5  =  0). 
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As  time  progresses,  the  front  moves  downstream  into  a  region  with  a  non-zero  vorticity 
disturbance.  Therefore,  the  numerical  solutions  include  transient  behavior.  The  relaxation 
to  steady  state  depends  on  both  the  acoustic  and  the  vorticity  responses  as  well  as  the 
angle  of  incidence.  If  9  is  above  the  critical  angle,  then  the  predicted  acoustic  response 
is  constant;  therefore,  steady  state  is  indicated  by  a  broad  flat  response  behind  the  front. 
If  9  is  below  the  critical  angle,  an  exponentially  decaying  acoustic  response  is  predicted; 
therefore,  steady  state  is  indicated  when  the  value  of  the  acoustic  response  at  the  front 
remains  essentially  constant  as  the  front  propagates  a  sufficient  distance.  Since  the  vorticity 
response  is  constant  in  both  cases,  steady  state  is  indicated  when  a  broad  flat  vorticity 
response  behind  the  shock  has  been  achieved. 

For  many  of  the  calculations  involving  either  large  disturbance  amplitudes  or  angles 
of  incidence  near  critical,  the  solution  did  not  asymptote  to  a  steady  state  even  after  the 
detonation  front  had  penetrated  a  significant  distance  into  the  unburnt  mixture.  The 
nature  of  this  transient  behavior  is  important  since  real  turbulence  consists  of  sudden 
non-steady  phenomena,  not  steady  plane  waves.  In  general,  the  transient  depends  on 
the  suddenness  in  which  the  front  meets  with  the  full  disturbance,  the  amplitude  of  the 
disturbance,  and  the  angle  of  incidence.  To  resolve  the  portion  of  the  transient  caused  by 
the  sudden  encounter  of  the  front  with  a  non-zero  disturbance,  the  calculations  allow  for 
a  slow  or  fast  ramping  of  the  amplitude  of  the  imposed  disturbance;  the  amplitude  as  a 
function  of  time  is  given  by 


€{t)  = 


6o(Mot/i?)2(3  -  2Mot/R) 
eo 


t  <  R/Mo 
t  >  RfMo- 


(19) 


Therefore,  when  the  position  of  the  front  is  x,  =  R,  the  flow  ahead  of  the  front  will  have 
its  maximum  disturbance  amplitude. 

The  value  of  the  exothermicity  parameter  a  is  chosen  to  be  equal  to  two  for  all 
calculations.  This  value  is  chosen  since  linear  theory  predicts  that  the  value  of  5  will  be 
significantly  different  from  its  value  in  the  non-reacting  case  for  almost  all  values  of  angle 
9  and  the  value  of  P  will  be  significantly  different  for  values  of  the  angle  near  critical. 
The  strength  of  the  detonation  is  chosen  so  that  the  normal  Mach  number  of  the  flow  is 
1.5  times  the  Chapman-Jouget  number.  The  critical  angle  for  these  parameter  values  is 
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=  24.89 


The  pressure  and  vorticity  responses  to  a  1%  disturbance  at  an  angle  of  ^  =  40°  for 
ramping  parameter  R  =  6.0  and  R  —  0.5  are  shown  in  Figures  3a  and  3b,  respectively.  For 
these  calculations,  32  grid  points  in  the  y-direction  are  used.  For  slow  ramping,  the  gas 
behind  the  front  responds  to  each  small  increment  in  the  amplitude  of  the  disturbance, 
and  a  smooth  approach  to  steady  state  is  observed.  For  fast  ramping,  R  =  .5,  the  acoustic 
response  shows  significant  overshoot  of  the  predicted  linear  value.  The  vorticity  response  is 
somewhat  smoother.  The  influence  of  heat  release  on  the  transient  is  considered  by  letting 
a  =  0,  R  =  .5  (Figure  3c)  and  prescribing  the  same  normal  Mach  number;  the  overshoot 
and  fluctuations  of  the  acoustic  response  are  considerably  less  than  for  the  reacting  shock, 
but  the  vorticity  response  shows  essentially  the  same  behavior  as  in  the  detonation. 

As  seen  from  Figure  3a,  when  the  ramping  is  slow  and  the  amplitude  of  the  disturbance 
is  small,  the  predictions  of  linear  theory  agree  with  the  calculated  acoustic  response  but 
there  is  a  discrepancy  in  the  vorticity  response  of  about  2%.  As  mentioned  in  Zang, 
Hussaini  k  Bushnell^’*^,  since  the  vorticity  response  involves  computation  of  derivatives  in 
both  the  X  and  y  directions,  the  vorticity  response  calculations  are  possibly  less  accurate 
than  the  pressure  response  calculations.  When  only  16  grid  points  are  used  in  the  y 
direction,  the  discrepancy  in  the  vorticity  response  increases  to  4-5%  and  the  pressure 
response  still  agrees  with  the  linear  predictions.  It  should  be  noted  that  the  calculated 
vorticity  response  is  always  less  than  the  predicted  values.  Another  important  parameter 
which  influences  the  transient  behavior  is  the  amplitude  of  the  disturbance.  For  the  10% 
disturbance  at  9  =  40°  (Figure  3d),  the  .'i.:oustic  response  takes  a  longer  time  to  asymptote 
to  a  steady  state. 

The  comparison  between  linear  theory  and  numerical  simulations  for  other  angles  is 
seen  in  Figure  4.  The  circles  are  the  calculated  responses  to  a  1%  disturbance  and  the 
diamonds  are  the  calculated  responses  to  a  10%  disturbance.  As  previously  found  for 
the  interaction  of  a  non-reacting  shock  with  a  vorticity  disturbance,  the  linear  predictions 
and  the  calculated  acoustic  responses  agree  for  disturbance  amplitudes  up  to  10%  when 
the  angle  of  the  incoming  flow  is  not  near  the  critical  angle.  The  calculated  vorticity 
responses  are  consistently  about  2%  below  the  predicted  values  when  32  grid  points  are 
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used.  Although  many  calculations  were  made  for  small  amplitude  disturbances  (1%  or 
less)  and  angles  within  5°  of  critical,  the  results  are  not  given  in  Figure  4  since  a  steady 
state  value  could  not  be  reliably  determined.  The  results  of  these  calculations  are  discussed 
below. 

For  disturbance  amplitudes  larger  than  10%,  the  solution  did  not  reach  a  steady 
state.  A  typical  run  is  shown  in  Figure  5  for  a  30%  disturbance  at  0  =  40°  with  ramping 
parameter  R  —  3.  Linear  theory  is  a  good  predictor  of  the  scale  for  the  overall  response, 
but  the  distortions  produced  in  the  detonation  front  prevent  the  solution  from  reaching 
steady  state;  in  fact,  the  distortions  prevented  the  calculations  from  proceeding  further 
since  the  time  step  based  on  the  CFL  number  became  exceedingly  small.  Furthermore, 
since  the  front  is  no  longer  planar,  the  comparison  between  the  least  squares  fit  to  the 
forms  (17)  and  the  linear  case  no  longer  holds.  The  same  behavior  was  observed  for  10% 
disturbances  at  angles  within  5°  of  critical. 

Linear  theory  predicts  that  exothennicity  has  its  greatest  effect  on  both  the  vorticity 
and  pressure  responses  near  the  critical  angle;  therefore,  it  is  important  to  resolve  the 
behavior  of  the  solution  for  angles  near  critical.  .4s  the  angle  of  incidence  approaches  the 
critical  angle,  it  is  observed  that  the  relaxation  of  the  transients  takes  longer  and  longer. 
These  transients  are  not  related  to  the  ramping  of  the  disturbance  amplitude  as  was  shown 
for  d  =  40°.  The  response  to  a  1%  disturbance  with  d  =  30°  is  shown  in  Figure  6  (16 
grid  points  in  the  j/-direction  and  ramping  parameter  i?  =  6).  The  transient  nature  of  the 
solution  provides  for  considerable  overshoot  of  the  predicted  acoustic  response  with  very 
little  relaxation  over  the  duration  of  the  run.  To  continue  the  calculations  past  .r,  =  10 
would  require  adding  additional  grid  points  in  the  x  direction  to  keep  the  same  resolution 
making  each  time  step  considerably  more  expensive.  In  addition  to  the  large  pressure 
fluctuations,  the  vorticity  response  also  displays  some  overshoot  which  is  not  observed 
for  slowly  ramped  disturbances  at  angles  more  than  5°  from  critical.  This  is  particularly 
surprising  since  calculations  at  other  angles  indicate  that  the  vorticity  response  remains  at 
5%  below  the  predicted  value  and  it  is  quite  smooth  for  such  slowly  ramped  disturbances. 
The  same  parameter  values  were  used  but  with  fewer  grid  points  in  the  x-direction  allowing 
integration  to  continue  until  x,  =  24.  Although  the  resolution  is  less,  a  relaxation  of  the 
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transients  is  observed;  the  acoustic  response  oscillates  about  a  value  accurately  predicted 
by  linear  theory,  and  the  vorticity  response  oscillates  about  a  value  which  is  approximately 
5%  below  the  value  predicted  by  linear  theory;  there  appears  to  be  a  slow  decrease  in  the 
amplitude  of  these  fluctuations.  These  results  are  consistent  with  steady  states  at  larger 
angles  of  attack  found  using  16  grid  points.  A  similar  trend  in  the  vorticity  response  for 
angles  of  attack  slightly  less  than  critical  develops.  The  vorticity  response  displays  a  slight 
overshoot  of  the  value  predicted  by  linear  theory  which  is  not  present  for  angles  away 
from  critical,  and  returns  to  a  small  oscillation  about  a  value  below  the  linearly  predicted 
value.  The  overshoot  of  the  predicted  pressure  and  vorticity  respoiuses  is  observed  for 
angles  of  attack  between  two  and  five  degrees  of  the  critical  angle;  however,  relaxation  of 
the  transients  is  significantly  slower  as  the  angle  of  incidence  nears  critical.  For  6  =  24.9° 
(Figure  7).  it  can  b^  ,een  that  the  acoustic  response  does  not  reach  the  value  predicted  by 
linear  theory  but  is  still  increasing  when  the  calculation  is  stopped.  The  vorticity  response 
is  also  still  increasing,  having  reached  only  half  of  the  value  predicted  by  linear  theory. 
This  calculation  was  terminated  due  to  the  increasing  expense  of  maintaining  an  accurate 
numerical  resolution;  however,  the  calculations  were  repeated  using  a  very  coarse  grid, 
and  it  was  found  that  even  though  the  pressure  and  vorticity  responses  continue  the  slow 
increase,  the  pressure  has  achieved  only  40%  of  the  linear  theory  value  and  the  vorticity 
only  60%  of  the  linear  theory  value  when  x,  =  50.  These  results  are  only  qualitative  at 
best  but  seem  to  indicate  that  the  lengthening  of  the  transient  response  for  angles  within 
2°  of  critical  counteracts  the  transient  overshoot  of  the  pressure  and  vorticity  responses 
found  for  other  angles  within  5°  of  critical. 

IV  CONCLUSIONS 

Nonlinear  calculations  of  the  response  of  an  initially  plane  detonation  wave  to  a  vor¬ 
ticity  disturbance  show  that  the  results  of  steady-state  linear  theory  are  useful  in  providing 
an  overall  scale  of  the  response.  In  cases  in  which  the  angle  of  incidence  is  near  critical, 
disturbance  amplitudes  are  moderate,  a:; d/or  there  is  a  sudden  encounter  with  a  distur¬ 
bance,  the  calculated  responses  display  a  transient  overshoot  of  the  linear  prediction.  It  is 
found  that  exothermicity  increases  the  overshoot  in  addition  to  increasing  the  value  of  the 
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predicted  linear  response.  Also,  the  significant  departures  of  the  predicted  responses  of 
the  reacting  shock  from  the  non-reacting  shock  near  the  critical  angles  appear  to  be  real, 
and  the  calculated  responses  show  transient  overshoot  of  the  predicted  values  for  angles 
between  two  and  five  degrees  of  critical.  Closer  to  the  critical  angle,  the  lengthening  of  the 
relaxation  time  for  the  transient  produces  a  competing  effect.  Also,  it  is  found  that  criti¬ 
cal  angle  of  linear  theory  is  an  accurate  predictor  of  the  transition  in  the  behavior  of  the 
acoustic  response.  Previously^’'*,  it  was  reported  that  the  change  from  a  constant  pressure 
response  to  an  exponentially  decaying  pressure  response  occurred  at  an  angle  significantly 
different  than  predicted  by  linear  theory.  For  the  situation  presented  here,  it  is  determined 
that  the  long  relaxation  of  the  transients  for  angles  slightly  above  critical  makes  the  pres¬ 
sure  response  appear  to  be  exponentially  decaying;  when  in  fact,  if  the  calculations  are 
of  a  sufficient  duration,  the  pressure  response  eventually  changes  and  eventually  oscillates 
about  a  a  constant  value  which  is  consistent  with  the  value  predicted  by  linear  theory. 

Since  exothermicity  is  seen  to  increase  both  the  vorticity  and  acoustic  responses  for 
the  detonation  wave/vorticity  wave  interaction,  it  would  be  appropriate  to  study  the  in¬ 
teraction  of  a  detonation  wave  with  a  fully  developed  vortex.  Such  a  numerical  study  was 
presented  in  Meadows,  Kumar  and  Hussaini®  for  the  non-reacting  shock  using  a  shock 
capturing  scheme  which  is  more  appropriate  than  the  shock  fitting  scheme  for  the  study  of 
shock  interaction  with  large  disturbances.  Of  particular  interest  is  the  ability  to  capture 
secondary  shocks  which  agree  with  experimentally  observed  features  of  the  flow.  This  is 
not  possible  using  shock-fitting  methods. 
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APPENDIX; 

The  improved  shock  fitting  scheme  utilized  by  Kopriva,  Zang,  and  Hussaini^  is  gener¬ 
alized  to  allow  for  exothermicity  of  the  reacting  shock  and  to  allow  for  a  time  dependent 
flow  eihead  of  the  front.  The  generalized  Rankine-Hugoniot  conditions  for  a  stationary 
reacting  shock  can  be  simply  written  as 


P2  =  Pi  +  g(M) 

(20) 

m  =  f{M) 

(21) 

g{M)  =  In - 

rY> 

(22) 

m 


where  the  subscript  1  refers  to  known  quantities  ahead  of  the  front  and  subscript  2  refers 

to  the  corresponding  quantity  behind  the  front;  f{M)  is  given  by  (8).  M  is  the  normal 

Mach  number  of  the  flow  ahead  of  the  front,  and  m  is  the  normal  velocity  behind  the  front 
divided  by  the  sound  speed  ahead  of  the  front.  If  all  other  velocities  are  nondimensionalized 
with  respect  to  the  speed  of  sound  for  the  mean  flow  ahead  of  the  front,  then  the  prescribed 
flow  has  a  local  speed  of  sound  ci  at  the  front  which  differs  from  unity  and  depends  upon 
both  y  and  t.  Denoting  the  front  position  by  fj  (y,  t)  ,  the  speed  of  the  front  traveling  in 
the  positive  x  direction  is  ^rj{y,t)  —  (17;, 0).  The  velocities  of  the  gas  in  the  frame  of 
reference  where  the  front  moves  Qi  =  are  related  to  M  and  m  by 

ciM  =  UjN,  -Ql  -N  (23) 

cim  =  UfN^  -Q2  - N  (24) 

where  N  =  (TV, ,  TV,  )  is  the  normal  to  the  shock  front  pointing  in  the  direction  of  the  domain 
subscripted  by  1.  The  key  to  the  shock  fitting  routine  is  the  compatibility  equation, 

P2t  -t-  -T^  •  Q2.  =  -Q  •  VP  -  -T^  •  R  -  cT^  •  VP  -  Rp  (25) 

c  c 

where 

R  =  (uu,  +  vuy ,  uvg  -I-  vvy)  Rp  =  7u,  -f  71;,  (26) 

derived  from  the  Euler  equations.  By  differentiating  the  Rankine-Hugoniot  relations  and 
equations  (23)and  (24),  an  equation  for  the  acceleration  of  the  detonation  front  can  be 
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found  and  is  given  by 


Ui,{y,t)  =  {-cicPu  +  cicC +  {cG--fF){Q^,-N +  Qi-N,) 

-(cG  -  7F  +  7Ci)f//iVj  t  +  C17Q2  • (27) 

+cit  [(cG  —  ')F)M  +  ci7m]}  {(cG  —  7F  +  '^ci)Nx  }”  ^ 

where  G  =  g'{M),  F  =  f'{M),  and  C  is  the  right  hand  side  of  equation  (25).  During 
the  calculations,  the  quantity  C  is  evaluated  using  the  solution  of  the  Euler  equations 
from  the  previous  time  step,  and  all  other  quantities  are  calculated  using  the  appropriate 
jump  conditions.  The  front  position  and  velocity  are  updated  for  each  time  step  using  the 
two-step  Lax-Wendroff  scheme. 
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Figure  1.  Schematic  of  model  problem  in  the  physical  domain.  Velocity  of  unbnnit 


mixture  is  IVo  + 
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Figure  2.  Schematic  of  model  pioblein  for  linear  theory.  Acoustic  disturbances  in 
the  burnt  gas  propagate  at  angle  to  the  x-axis.  Vorticity  and  Entropy  disturbances  are 
convected  with  tlie  mean  flow  at  angle 
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Figure  3a.  Presiure  and  Vorticity  teaponse  to  a  1%  tlowly  tamped  disturbance. 
A  =  6,  a  =  2,  and  0  =  40*.  Solid  tine  is  linear  theory.  The  curves  correspond  to 
t  =  .380,  .613,  .847,  1.140,  1.341. 
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